A general T-matrix approach applied to two-body and three-body problems 

in cold atomic gases 
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We propose a systematic T-matrix approach to solve few-body problems with s-wave contact 
interactions in ultracold atomic gases. The problem is generally reduced to a matrix equation 
expanded by a set of orthogonal molecular states, describing external center-of-mass motions of 
pairs of interacting particles; while each matrix element is guaranteed to be finite by a proper 
renormalization for internal relative motions. This approach is able to incorporate various scattering 
problems and the calculations of related physical quantities in a single framework, and also provides 
a physically transparent way to understand the mechanism of resonance scattering. For applications, 
we study two-body effective scattering in 2D-3D mixed dimensions, where the resonance position 
and width are determined with high precision from only a few number of matrix elements. We also 
study three fermions in a (rotating) harmonic trap, where exotic scattering properties in terms of 
mass ratios and angular momenta are uniquely identified in the framework of T-matrix. 
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I. INTRODUCTION 

Interacting ultracold atoms have gained a lot of re- 
search interests for their interaction strength and dimen- 
sionality are highly controllable by making use of Fesh- 
bach resonance and external confinements jl], Q ■ In such 
dilute atomic gases, the interaction between atoms can be 
well approximated as contact potential which is charac- 
terized by the s-wave scattering lengthp], In this con- 
text, the few-body problems play very important roles 
in studying many-body properties. For instance, solu- 
tions of these problems determine effective interactions 
between atom-atom, atom-dimer and dimer-dimer[3, 0], 
which are fundamental elements to formulate the many- 
body effective Hamiltonian; moreover, the consideration 
of two-body short-range physics leads to a series of exact 
universal relations for a many-body system, as first pro- 
posed by Tan0-0| and recently verified in experiment [§[■ 

Previous studies of few-body systems have revealed 
many nontrivial effects. One famous example is the Efi- 
mov effect for three atoms depending closely on 

the short-range interacting parameter apart from sin- 
gle s-wave scattering length. Another typical effect is 
the two-body induced resonance scattering and induced 
bound state due to external confinements [Tl]-r24j . Among 
most of previous studies, the problems were solved in 
the framework of two-channel models 25| or by using 
Pseudopotentials[2a. |27|. In this article, we present us- 
ing T-matrix approach to solve few-body problems with 
s-wave contact interaction in the field of ultracold atoms. 
Compared with other methods, T-matrix is able to sys- 
tematically provide exact solutions for few-body prob- 
lems, and more importantly, is able to work in a much 
efficient and physically transparent way. 

In this article, we shall first formulate T-matrix 
method and introduce its essential concept, i.e., the 
renormalization idea to integrate out all high-energy(or 
short-range) contributions for relative motions. Then a 
systematic treatment is presented to a general N-body 
system with contact interactions and with possibly trap- 



ping potentials. The key point of this method is to make 
use of the interaction property and introduce a set of 
orthogonal molecular states, which describe the external 
center-of-mass(CM) motions of pairs of interacting par- 
ticles; then the problem is generally reduced to a matrix 
equation, and a proper renormalization scheme for the in- 
ternal relative motions ensures finite value of each matrix 
element. Using this method, we obtain the bound state 
solution, scattering matrix element and reduced coupling 
constant in the low-dimensional subspace. Moreover, the 
treatment has a lot more physical meanings and allows 
us to make analytical predictions to various induced res- 
onances under confinement potentials. In all, T-matrix 
approach is able to unify many studies of different issues 
in the single framework. 

To show the efficiency of this approach, we apply it 
to study two-body effective scattering in 2D-3D mixed 
dimensions, where T-matrix not only provides a trans- 
parent way to understand the mechanism of multiple 
resonances, but also gives explicit expressions for the 
resonance position and the width. Particularly, its ef- 
ficiency lies in that each resonance can be determined 
accurately by only calculating a few number of matrix 
elements. Moreover, we apply this method to study three 
two-component fermions in a (rotating) harmonic trap, 
and show its unique advantage in identifying scatter- 
ing properties in different angular momentum channels. 
The ground state is obtained for a rotating and trapped 
system, which gives important hints for quantum Hall 
physics in the fermionic atom-dimer system. 

The rest of this paper is organized as follows. In sec- 
tion II, we introduce the renormalization concept and 
present T-matrix formulism to solve a general few-body 
problem. Its relations to other approaches, advantages 
and limitations are also discussed. Section III is the ap- 
plication of T-matrix approach to two-body problems, 
where specifically we study the scattering resonances in 
2D-3D mixed dimensions. Section IV is the application 
to three fermions in a (rotating) harmonic trap, where 
the energy spectrum and scattering property are studied 
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for different mass ratios and angular momenta of three 
fermions. We summarize the paper in the last section. 



II. T-MATRIX APPROACH 

In this section we give a systematic prescription of T- 
matrix approach to solve few-body problems. The re- 
sulted matrix equation is given by EqJTTJin Section IIA, 
from which we extract three important physical quanti- 
ties as given by Eqs. (|12ll4ll8j) in Section IIB. We also 
discuss its relation to other widely used methods in Sec- 
tion IIC, and demonstrate its unique advantages and lim- 
itations in Section IID. 



A. Basic concept and general formulism 

We start from the Lippmann-Schwingcr equation based 
on standard scattering theory, 



= \iJo) + Go(E)T\iPo). 



(1) 



Here \ipo) is the eigen-state for non-interacting system, 
with Hamiltonian Hq = + Vr composed by kinetic 
term and external trapping potential; is the scattered 
state in the presence of interaction potential U; 
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H Q + iS 



(2) 



is the Green function; the scattering matrix can be ex- 
panded as series T = U + UG U + UG UG U + ... which 
leads to 



T 



(l-UGo^U. 



(3) 



To this end, Go, T and U are all matrixes expanded by 
certain complete set of states. If we use {(V'o)} to ex- 
pand T-matrix, then each T-matrix element in EqJ5] di- 
rectly relates to the scattering amplitude in the scattered 
wavefunction and therefore represents the effective inter- 
action in the low-energy space. To obtain the effective 
interaction, a physically insightful way is to employ the 
concept of renormalization. For the fundamental two- 
body s-wave scattering with contact interaction, 



U(r) = U Q S 3 (r), 



(4) 



one can resort to a simple momentum-shell renormaliza- 
tion scheme 20] . The spirit here is to take all intermediate 
scattering from low-k(momentum) space to the shell in 
high-k space as virtual processes, which in turn modify 
the effective interaction strength(CZ) in the low-k space 
by a small 5U; perturbatively 5U in terms of the shell 
momentum((5A) follows 



A-<5A<|k|<A 



(5) 



here V is the volume, and £k = k 2 /(2^) is the energy 
for relative motion of two particles with reduced mass fi. 
The resultant RG flow equation reads 28] 
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(6) 



relates the bare interaction Uq to the zero-energy effective 
one(Xo = 27ra s //i) via, 



27ra, 



1 1 V 1 
U~ + V V^k' 



(7) 



with o s the s-wave scattering length. For a general N(> 
2)-body problem, the idea of renormalization, though not 
as explicitly shown as above, always serves as the under- 
lying principle through the whole scheme(see below). 

Now we proceed with the general T-matrix approach. 
Suppose a Q— species system, and the i-th (i = 1...Q) 
species has Ni identical particles residing at x^,...x^. ; 
Ui and Uij{i < j < Q) are respectively the bare interac- 
tion strengths between particles within the i— th species 
and between different species (z and j); Ui (£/y) is re- 
lated to the corresponding scattering length at (ay) via 
Eq0 Taking advantage of the zero-range property of the 
interaction 



Q n z 



i—l m<n 



i<j m—1 n—1 



and thus the same property of T-matrix given by Eq|31 we 
expand U and T by a set of molecular states { |x^ — yJ n — 
0,A)}. This state is defined such that one pair of in- 
teracting particles (x'„, x£) locate at the same site, and 
A is the energy-level index for the remanent degrees 
of freedom. More detailed description of the molecu- 
lar state is given in Appendix [A] Further, for identi- 
cal bosons/fermions the molecular state should further 
be symmetrized/antisymmetrized by superpositions of 
above individual ones. Explicitly we have 



^o)=^/ A J |r-/ = 0,A) ) 



(9) 



IX 



here each state, |n = 0, A)(rj = |x* n — x£|) with I < 
max(J) = ; corresponds to one interaction term 

in EqJSJ The coefficients {//} satisfy 



i'\ r 



hi 
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Sii>S xx , - C{{,\ = (n = 0, A|Vo>(10) 



with 



nil 

°AA< 



= ^Y,—6n'6xx> + (rj = 0,X\G \rr = 0,A'). 



(11) 

Detailed derivation of Eq[TU] is given in Appendix [Bj 

For the two-particle scattering in free space, the en- 
ergy level {A} is characterized by the CM momentum, 
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which is conserved by the interaction and thus irrele- 
vant to the scattering problem for relative motions. The 
molecular state is then as simple as \r = 0), with r 
the distance between two particles. EqO is then re- 
duced to T\ip ) — f\r — 0), and / is given by EalTOl 
as f^ 1 oc aj 1 + ik, reproducing the well-known relation 
between the scattering amplitude and the s-wave scat- 
tering length (<z s ). The applications of this approach to 
other few-body systems will be introduced in Section III 
and IV. 



B. Calculation of physical observables 

With the information of molecular states, all physical 
quantities can be deduced straightforwardly. We shall 
enumerate below three quantities that are detectable or 
observable in experiments. 

(I)Bound state solution. In this case |^»o) is absent, 
Eqfiniis given by the pole of T-matrix(see EqJ3]), i.e., 



Det(l - UG (E b )) = 0, 



(12) 



where E b is the binding energy, and the eigen- vector {/^} 
gives the bound state as 



(13) 



IX 



For two-particle scattering(with scattering length a s ) in 
free space, Eqs. (fT2l H3l) give E b = —1/(2)10%) and \ipb) oc 

E 



(II)T-matrix element. Generally, T-matrix element 
between \i/jo) and itself characterizes the scattering prop- 
erty of low-energy particles. With Eqs. (|9ll0[) . we obtain 



(14) 



IX 



which can also be obtained from Eq[3] as ("001(1 — 
UGo)~ 1 U\ipo). For two-particle scattering in free space, 
fl = fis proportional to the scattering amplitude. 

(Ill)reduced interaction. If trapping potentials con- 
fine atoms in a lower dimension, there are two distinct 
scattering channels for the low-energy state, namely the 
open(P) or closed(Q) channel, depending on whether its 
wavefunction propagates or decays at large interparti- 
cle distance in the lower dimension. The effective in- 
teraction strength for low-energy particles in the open 
channel is modified from the original bare one by vir- 
tual scatterings to the closed channels. We assume g e g 
as the modified interaction strength in open channel, 
which, for instance, has been defined in the reduced ID 
Hamiltonian [12], ll3| |2lj, [29( under tight transverse har- 
monic traps. Following the same procedure in obtaining 
T-matrix element in (II), we have 



9cS = Upp + UpqGqUqp + UpqGqUqqGqUqp 



(15) 



where Ut s (t,s — {P, Q}) is the bare coupling strength 
between two particular channels; 



G, 



Q 



1 



E-H, 



Q 



id 



(16) 



is the Green function for Hamiltonian Hq that only acts 
on closed-channel states. Again due to the property of 
[/(EqJS]), we insert into Eq[T5] a set of molecular states 
Eja \ r i = 0,A)(rr = 0,A|[30j. Assuming two column 
vectors 

£ = {(r 7 = O,A|0 o )}, C={v-z = O,A|*70o)}, (17) 



we then obtain 



9cs 



ec+e(uG%)c+exuG$) 2 t+ 



ea - uc$y 

dj \(l-UG$ 



)- 1 C/|^ > 



(18) 



here UGq is a matrix expanded by molecular states. 
Note that g e g is different from (^ol^l^o) in (II) only by 
the Green function therein. Obviously g e g is the renor- 
malized coupling strength in the open channel by all vir- 
tual scattering to states in closed channels. 



C. Relations with other methods 

In this section, we analyze the intrinsic relation be- 
tween T-matrix and other widely used methods, such as 
those in the framework of two-channel models [13, [lj| and 
pseudopotentials [i~i"14i~3l Il8l . 

In two-channel models, the closed-channel molecules 
are explicitly included in the Hamiltonian; these 
molecules couple to atoms in open-channel and thus me- 
diate interactions between the atoms. In the present T- 
matrix method, the molecular states(defined in Eq|H]) can 
be considered as the analog of closed-channel molecules 
in two-channel models. The similarities lie in that they 
are both constructed in a way that follows the zero-range 
property of interaction, and describe the CM motion of 
two interacting particles. 

In pseudopotentials, the problem is solved by applying 
Bethe-Peierls boundary condition to the wavefunction at 
short inter-particle distance, i.e., 



lim 0(xi,x 2 ,x 3 ...) = (- -)/(xi 

r / = |xi-x 2 |-vO r/ a I 



x 2 ,x 3 ...), 
(19) 

here the index / denotes the pair {xi, x 2 } and aj is the 
scattering length between them. On the other hand, we 
notice that above asymptotic behavior can be automati- 
cally satisfied by the present scheme of T-matrix method. 
To show this, we examine EqfT]at short inter-particle dis- 
tance by projecting it to certain molecular state, 

lim (r7,A|V-) - (77 = 0, A|0o) + 
n— s-o 
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J2fy{ri = 0,MGo\rr=0,X) 

k 

= lim -I). (20) 

To derive Eq|20] we have used EqJTO] and the Fourier 
transformation of zero-energy Green function in free 
space. Remarkably, Eqj20] shows that the asymptotic 
form is given by the inverse of bare potential 1/Ui, which 
is universal regardless of any trapping potential. In this 
sense, the Bethe-Peierls boundary condition fEq fK))) in 
the framework of pseudopotentials is equivalent to renor- 
malization equation (Eq(7]) in T-matrix method. 

D. Advantages and Limitations 



ordinate method in unitary limit [H, H3| , or by emplo yin g 
a three-body force to eliminate the cutoff dependence[35| . 
The extension of the present T-matrix approach to iden- 
tify such nontrivial few-body effects is out of the scope 
of this paper. 



III. APPLICATION TO TWO-BODY PROBLEM 

In this section we consider the two-body system. First 
we present the formulism of physical quantities intro- 
duced in Section IIB, for both cases when the CM motion 
and relative motion can or cannot be decoupled. Finally 
we apply this method to address the effective scattering 
of two particles in 2D-3D mixed dimensions. Particularly 
we emphasize the physical insight given by T-matrix to 
understand the resonance mechanism, as well as its effi- 
ciency in realistic calculation of scattering parameters. 



In principle, the T-matrix fomulism presented in Sec- 
tion IIA is appliable to a general few-body problem with 
zero-range interactions. It gives a unified treatment to 
different scattering issues, as shown in Section IIB. Be- 
sides, T-matrix approach has the following advantages. 

First, we compare this scheme, using molecular states 
to expand T-matrix fEq llOl with that using original 
N— particle states. The latter requires a matrix dimen- 
sion as large as T N (T is the cutoff index for single par- 
ticle energy levels), while the former could reduce it to 



at most 
Q(Q+i) r N 

2 



Q(Q+i) r N-i 



for a general case and further to 
for the special case when CM motion can 
be separated out. 

Second, this scheme is physically insightful in that it 
reveals general couplings between the CM and relative 
motions. Specifically, CM serves as external indices of the 
matrix, while relative motions contribute to each element 
by the renormalization of internal degree of freedom. As 
we shall see in Section III, this picture is essentially im- 
portant to understand the mechanism of resonance scat- 
tering in the two-body system. 

More advantages related to the realistic calculations 
will be shown in Section IIIB and IV, when applying this 
method to specific two-body and three-body problems. 

However, the present T-matrix method still has lim- 
itations under certain circumstances. When examining 
Eqs. (|10llip , one can see that the convergence of the so- 
lution from the matrix equation generally requires two 
conditions. First, each matrix element is finite; secondly, 
the solution is independent of the matrix size. Any viola- 
tion of above two conditions implies that the zero-range 
model is insufficient to characterize the interacting sys- 
tem, such as when Efimov physics emerges in the three- 
body sector d, [13. T-matrix is able to identify the vio- 
lation of the first condition (see Section IV). However, it 
would be quite involved for it to identify the second con- 
dition. Alternatively, when Efimov physics appear, one 
can solve the problem by resorting to hyperspherical co- 



A. Formulism 

1. r and R, decoupled system 

For trapping potential Vr{xi,x 2 ) = Vr(R) + Vr(r), 
the molecular state is simply \r = 0, Ao), with Ao char- 
acterizing the CM motion and staying unaffected by the 
interaction. The bound state solution Ej, is determined 
from a single equation as 



/' 



27ra< 



- C(E b ), 



with 



C(E) 



E-Ei+iS' 



(21) 



(22) 



here {4> n (r)} is a complete set of eigen-states only for 
the relative motion. A convergent solution of Ef> requires 
that the ultraviolet divergence in each term of Eq[22] be 
exactly cancelled with each other. This is actually satis- 
fied by a regular potential Vr(x) without singularity at 
any position x(36[. 

The T-matrix element, T mn = (m\T(E)\n) , which rep- 
resents the scattering amplitude from initial state \n) to 
final state |m), is given by 



</4(o)0„(o) _ [i 



-cm- 



(23) 



Compared with Eq|2TJ this shows that a bound state 
emerges when all elements T mn simultaneously diverge. 

Similarly for reduced interaction in the lower dimen- 
sion, we obtain 



\MQ)\ 2 

9cB 



/' 



27ra s 



C Q (E), 



(24) 
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where follows the form of Eqj22] with the second 
summation over all closed channel states. In this sense 
the confinement induced resonance(CIR), referring to 
<?eff 00 7 occurs at 



According to T- matrix method, assuming Eg iTTI yields 



C Q (E). 



(25) 



This in turn determines a closed channel bound state 
with the same energy as \ipo). Compared to existin g ex - 
ploration of CIR mechanism in quasi- ID system [T2T 1 13j. 
T-matrix shows in a general way how the divergence of 
reduced interaction in the open channel is associated with 
the emergence of a closed-channel bound state at the 
same energy level. 

Finally, Eqs. (|22|23|24|) indicate the way how the trap- 
ping potentials modify the low-energy scattering the- 
ory. The modification is in fact through the intermedi- 
ate virtual scattering processes, i.e., by redistributing the 
energy levels and changing coupling strengths between 
these states. Above formula can be applied to ordinary 
harmonic confinements studied before |Tl| - [l3l [l5T - [l7] . 



2. r and R coupled system 

For a general trapping potential, the two-body non- 
interacting Hamiltonian can be divided to three pieces, 
describing the relative motion H re i(r), CM motion 
-ff cm (R), and couplings in-between H cp (r, R). The 
molecular state is then introduced as \r = 0,A), and 
$a(R) = (R|A) is the eigen-state of 

H cm {B) = -^ + V T>1 {Tl) + V T , 2 {n). (26) 

Combining with Eq llOl one can obtain all solutions cor- 
responding to (I, II, III) in Section IIB. 

In this case, the trapping potential can induce multiple 
two-body scatterin g re sonances as revealed previously in 
several settings [2ll |23|. l2~j ] by numerical calculations in 
coordinate space. Next we show that these resonances 
can be analytically figured out in the framework of T- 
matrix method. We classify the situations by whether 
the effective scattering is in 3D space or in the reduced 
lower dimension. 

When external trapping potentials are applied but the 
low-energy scattering wavefunction, ■)/ ; (xi,X2), can still 
behave in a propagating way at large 3D inter-particle 
separations |xi — x 2 | — > 00, then its asymptotic form can 
be written as 



0>eS 



d(xi,x 2 ) 



(27) 



Here d is the modified interparticle distance according 
to the confinement (see also Ref.[23| and discussions in 
Appendix [Cjl ; a c g is the effective scattering length, and 
can be directly related to T— matrix element for zero- 
energy scattering state [37 1. 



CaA'-^AV — Cv-F\vi 



(28) 



with c v (y = 1,2...) the eigenvalue and T v the corre- 
sponding eigenvector, then we have 



(M T m) = 2^ — m — 



(29) 



2iTa s 



This equation explicitly predicts an infinite number of 
resonances (a c g — > 00) when each discretized c„ individ- 
ually match with 2 £ a by tuning a s in realistic experi- 
ments. The resonance position and width can be con- 
veniently extracted from the exact diagonalization of C- 
matrix. 

To explore the mechanism of these resonances, first 
we only focus on the diagonal elements of C-matrix. 
Within each molecular channel, all relative motion levels 
are coupled together by the attractive interaction and 
this potentially leads to a bound state. This bound 
state(relative motion) combined with the molecular chan- 
nel(CM motion) tend to produce the zero total energy, 
and thus give rise to the divergent T-matrix or a e g — > 00. 
By tuning the interaction or a s , the zero-energy state 
will emerge in order from each molecular channel and 
cause the resonance of a c g. The width of each resonance 
is determined by the coupling between the zero-energy 
scattering state and each molecular state, which becomes 
narrower for higher levels of molecular states. 

However, above understanding of multiple resonances 
is not rigorous, because different molecular channels 
could also couple with each other by the combination 
of H re i(r) and H cp (r, R). This additional coupling, as 
shown by off-diagonal elements of C-matrix, would give 
a correction to the ideally predicted resonance position. 
In the following section, we shall address these issues by 
studying a specific system with multiple resonances, i.e., 
two atoms scattering in 2D-3D mixed dimension. 

When external trapping potentials are applied such 
that at low energies, V( x 1j x 2) only propagates as two 
particles are far apart in a lower dimension(open chan- 
nel), then the same analysis can be applied to the ef- 
fective interaction g e g in this channel. Now C-matrix is 
defined by the closed-channel Green function Gq (Eq|T5|) , 
which equally results in the matrix equation 



XX' jX' 



Q 

7 A- 



Combined with EqfT8l it gives 



^eff 



E 

7 



_B r Q 

2-na s L T 



(30) 



(31) 



Therefore g e g will go through a resonance as long as 
one is matched with 2 ^ a by tuning a s . This cor- 
responds to the energy of dressed bound state in each 
closed molecular channel moves downwards and touches 
the threshold energy of open channel. The resonance 
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width would be narrower for higher molecular channels 
due to the smaller overlap with the low-energy scattering 
state in open channel. 



B. Results of scattering in 2D-3D mixed 
dimensions 



We consider one atom ( 41 K or 40 K, labeled by A) is 
axially trapped by a tight harmonic potential with fre- 
quency lu a , while the other atom ( 87 Rb or 6 Li, labeled 
by B) is free in 3D space. The Hamiltonian reads 



1 



2m A '"2 W ^ 2m B 



V?B 1 U 6 3 (r A -r B ) 



H(r A ,r B ) 

mi \ l mi a 

(32) 

As shown in Appendix O the effective scattering length 
a e ff for zero-energy scattering is defined by the two-body 

wavefunction when d AB = \ ~^tPab + z b ~ * 00 ' 



with 



iI>(Pab>za,zb) <f>o(zA,a )(l - -7^-), (33) 

a-AB 



2na. 



V(il>o\T\i/>o) 



(34) 



Here 4> n (z,ao) is the eigen-state of ID harmonic os- 
cillator with characteristic length ao = \J\j (m A uj A ); 
/i = m A u/(l + u) is the reduced mass and u — m B /m A 
is mass ratio. 

The molecular state in this case is \r = 0,7V), where 
\N) denotes an eigen-state of the following Hamiltonian 



H cm (Z) 



1 



2(m A + m B ) dZ 2 2 



d 2 1 

+ -m A uj A Z 2 , (35) 



which describes the CM motion along (trapped) z di- 
rection, with oscillation frequency oJ — ui A j\J\ + u and 
characteristic length a — ao/(l + u) 1 / 4 . According to 
Eas. (|29l34p we obtain 



a c ff 
a 



5-/ lin 

N=0,2,... a s 



w N 



ejv 



(36) 



here ejv is the eigen- value of C— matrix fEq |C6l) deter- 
mined by Cnm^mn' = gn'Fnn 1 ] the resonance width 
is given by 



Wn — | ^ -^jV'jv/iV';0,o| 2 , 
N' 



(37) 



with f(N;n A ,k z ) defined by EqEU Due to the 
contact interaction and reflection symmetry of trap- 
ping potential, only molecular states with even parity 
(N = 0,2,4...) are relevant in this case. Fig. 1(a) and 
Fig. 2 (a) shows the first five resonances of a e s/ao for 
41 K- 87 Rb(rnA < m B ) and 40 K- 6 Li(m j4 > m B ) mix- 
tures, when tuning ao/a s from the weak coupling(— oo) 



to strong coupling(+oo) side. The + l)-th resonance 
of a e s is characterized by the position (ao/a s ) res = ejv 
and the width Wn- 

Amazingly, we find good accordance between e^v 
and each diagonal matrix elements Cnn, as shown by 
Fig. 1(b) and Fig. 2(b). That means the correction caused 
by off-diagonal couplings between different molecular 
channels are actually negligible. There are mainly two 
reasons for this. First, the amplitudes of off-diagonal 
C— matrix elements (Cnn') ar e much smaller than di- 
agonal ones, and decrease rapidly as \N — N'\ increases. 
Secondly, it can be attributed to the destructive inter- 
ference among couplings with different molecular chan- 
nels. To see this, we carry out a perturbative calculation 
in terms of off-diagonal couplings between neighboring 
molecular states, and get the relative correction to the 
N— th eigen-value as 



\Cnn< 



N'=N±2 



Cnn — Cn> 



-/a 



NN- 



(38) 



JV' 



Note that N — 2 and N + 2 terms contribute to 
with opposite signs, which results in further suppressed 
A^v- Insets of Fig. 1(b) and Fig. 2(b) show for the first 
eleven resonances; we see that the most significant effect 
of off-diagonal couplings occurs for the first resonance 
but is still negligibly small (~ 1.5% for K-Rb mixture 
and ~ 0.1% for K-Li mixture). 

The vanishing off-diagonal couplings between different 
molecular channels establish the unique advantage of T- 
matrix scheme, i.e., the resonance can be accurately de- 
termined by only a few number of related matrix ele- 
ments. In the limit of zero couplings, we have Tn'N ~ 
Sn>n and therefore ejv = Cnn, Wn = |/w ; o,o| 2 (see 
Eq |C4[) . In this limit and particularly for resonances at 
a s > side, the resonance position can be determined 
by matching bound state energy in each CM molecular 
channel with the threshold energy of two particles, i.e., 



(JV 



1, 



1 



UA 

2 ' 



(TV = 2,4,6...) (39) 



as shown by the dashed lines in Fig. 1(b) and Fig. 2(b). 
This equation, as is the direct outcome of T-matrix 
method, has been used previously to determine the reso- 
nance positions [Hj]. 

In addition, Fig.l and Fig. 2 show that the first 
resonance, which is mainly due to the coupling be- 
tween relative motion levels within the lowest molecular 
channel(JV = 0), always occurs at a s < side regardless 
of u > 1 or u < 1. The resonance position, however, 
sensitively depends on the value of u, as shown by the 
vertical dashed lines in Fig. 1(a) and Fig. 2(a). This phe- 
nomenon is closely related to the distinct behaviors of 
Cnn for different u. On one hand, when u <Si 1 we can 
omit the n^-dependence in In —function in Eq |C6[ then 
using Eq |C5l we obtain Cnn — > for any JV. This implies 
that when trapping the heavy atom, the resonance po- 
sitions [(ao/a s ) res ] tend to highly aggregate around un- 
shifted position(ao/a s = 0). On the other hand, when 
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(a) 





a/a. 



(a) 



-1.5 



-1.0 



-0.5 



0.0 
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0.5 



1.0 



1.5 




FIG. 1: (Color online) (a) .Effective scattering length a e ff/ao 
as functions of a /a s for 41 K(2D)- 87 Rb(3D) mixture. The 
dashed line denotes the first resonance at ao/a s = —0.43. 
(b)Diagonal matrix element Cnn of Eo lC6f denoted by x), 
the corresponding resonance position (ao/a s )res(+) and reso- 
nance width(o). Red dashed line is the function fit according 
to Eg 1391 Inset shows the relative correction A at as defined 
in Eg 1381 N = 0,2,4... in (b) respectively correspond to the 
(N/2 + I)— th induced resonance in (a) from left to right. 



u 3> 1, fN-nA.k- is vanishingly small for finite N, then 
the first term dominates in Eq |C61 This predicts res- 
onances approaching a s = CP. Only for large enough 
N the rest terms in Eq |C6l would dominate and predict 
resonances at a s > side. 

We note that the dependence of the first resonance po- 
sition (at a s < side) on the mass ratio u is in qualitative 
agreement with that for 0D-3D mixtures [2 2]. Actually 
we can gain the physical insight of such feature from the 
analysis of interaction potentials affected by the confine- 
ment. Suppose a square- well interaction potential U(r), 
which is — Vq at r < ro and zero otherwise, between A 
and B atoms. As Vb increases, the first scattering res- 
onance occurs at the critical value V c — (tt /2) 2 / '(2/irg) 
with n the reduced mass. When A or B is trapped and 
becomes localized, \i will be effectively enhanced, which 
reduces the critical V c and gives new resonance at a s < 
side, /j, and V c can be substantially modified if the lighter 
atom is trapped, and the resonance position will move 



1.5 




2 4 6 8 

N 



FIG. 2: (Color online) Same as Fig.l except for 40 K(2D)- 
6 Li(3D) mixture. The dashed line in (a) denotes the first 
resonance at ao/a s — —0.02. 

far away to a s = CP side. On the contrary, if the heav- 
ier atom is trapped, /x and V c would be little affected by 
the trapping potential, giving almost unshifted resonance 
near a s = oo. This analysis leads to similar conclusions 
for the resonance scattering in other mixed-dimensional 
systems, such as 1D-3D mixtures. 



IV. APPLICATION TO THREE-BODY 
PROBLEM 

Besides the two-body system, the general formulism of 
T-matrix approach allows its straightforward extension 
to other few-body systems. In this section we focus on a 
three-body system composed by two-component fcrmions 
in a (rotating) harmonic trap. We shall first present 
the formulism and then explore the interesting scattering 
property and identify the ground state level crossing in 
this system. 



A. Formulism 

We consider three fermions with one spin- J, (xi) and 
two identical spin-f (x2,xa) in an isotropic harmonic 
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trap. According to Eqs. (|A10HA13l) . we transform the here 
vector X = (^/2mjxx, -^/2m-j-X2, yj2m^^) to Y± = 

(V2M«R,V^r ± ,V^P ± ) by Y± T = A ± X T . Here 
(R, r_,p_) and (R, r + ,p + ) are all Jacobi coordinates, 



respectively corresponding to the effective mass Mr, /j,, fj,; 
the CM coordinate R and its mass Mr follow Eq lAllI 
the other Jacobi coordinates are 



x 2 - Xl, 



P 



with the same mass fi = 
reads 



1*3 



(40) 



the transfer matrix 




V ^/A/(rat+m ; ) V Af(m t +mj,) 



M 



J 

(41) 

we further obtain p + , r + by exchanging x 2 X3 in 
r_ , p , and obtain A+ by exchanging the second and 
third column of A- . 

Taking advantage of the property of transfer matrix 
(Eq lA14| ). we can see that with the same trapping fre- 
quency w, all three Jacobi coordinates (R, r± , p± ) can 
be well separated from each other. Independently one 
can also prove that the total angular momentum is 
also separable as Yli=i 2 3 F a (xi) — L a (R) + L a (p±) + 
£a( r ±) (a = x,y,z). Therefore for a trapped system 
with rotating frequency around z-direction, the rele- 
vant Hamiltonian in the rotating frame reads 



H(p ± ,r ± ) 
here 



H (p±) + H (r ± ) + U 6(r+) + U S(r_) 



H (r) = -^ + ^ 2 r 2 -QL z (r), 



(42) 



(43) 



The molecular state is defined with respect to Fermi 
statistics, 



= 0,A)-|r+=0,A)). 



(44) 



Here the first and second A represent the identical energy 
level {nlm} for the motions of p and p + under Hamil- 
tonian Ho. ({n} and {lm} are respectively the radial 
and azimuthal quantum number). Then we obtain the 
C— matrix element as 



C 



^ E - E x - E u + i5 



XX' 



with 



Fxx'i 
(45) 



-_ =0,A|G o |r+ = 0,A' 
r «(p)^ v (-/3p)~ 



^(O)^(-qp) 
E v + i 



^ E-E x 



(47) 



and a 2 + (3 2 — 1. To obtain EqHB] we have inserted 
into the Green function a complete set of eigen-states 
{v, A'} for the motions of (r + ,p + ). F; 



XX' 



Fy X here 



induce the coupling between different molecular levels, 
and non-zero F\x' require azimuthal quantum number 
{lm} be conserved. Note that the off-diagonal coupling of 
molecular states here is due to the many-body statistics, 
in contrary to the previous two-body case which is due to 
the external trapping potentials. More details regarding 
to the evaluation of Eqf45] are presented in Appendix [D] 



B. Results 

In the first part of this section we use T-matrix 
method to analyze the exotic scattering property of three 
fermions in different limits of mass ratios and in differ- 
ent angular momenta channels. In the second part, we 
present the energy spectrum and identify the energy level 
crossing between different angular momenta states for the 
(rotating) trapped system. 



1. Scattering property 



By analyzing Eqs. (|45l46p . we find nontrivial scattering 
properties at two limits of mass ratio u — m^/m^. Figj3] 
shows the schematic plots of Jacobi coordinates (r_ , p ) 
in both limits of u — > and u — > 00. 




FIG. 3: (Color online) Jacobi coordinates (r_,p_) for three 
fermions (fN*) in the limit of u — > (a) and u — > 00 (b). 
u = m-f/mj, is the mass ratio, a is given by EaliTl 



First, when u — > as shown by FigEJa), a — >■ 1, 
0, we have 

FxX' ~ SinSi' n. 



(48) 



Therefore the diagonal C-matrix for I ^ indicates the 
atom-dimer uncorrelated system with energy E — E a + 
Ed. Physically we can see from Fig(3Ja) that, the dimer 
formed by a light f and heavy J, is almost equivalent to 
single 4, atom, so the other f has s-wave interaction with 
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this dimer only when I — 0. Here the heavy I dominates 
the whole physics. 

Second, in the opposite limit when u oo as shown 
by FigCUb), the result is completely different. We find in 
this limit, 



C 



l^(0)| 2 



k 



E - E x - E v + i5 1 



-1)'])*AA' 

There are two direct consequences as follows. 

(i) for odd I, C\\> = oo, i.e., unphysical divergence in 
the high-energy space can not be properly removed. This 
is exactly the evidence of Efimov effect for large u where 
another short-range parameter is required to help fix the 
three-body problem [a. [Toj. 

(ii) for even I, U takes no effect and the sys- 
tem just behaves like non-interacting. This result is 
consistent with that obtained by Born-Oppenheimer 
approximation(BOA) [40j . Under BOA, the wavefunction 
is given by 

■0(xi,x 2 ,x 3 ) = kp(|x 2 -xi|) +7<^(|x 3 -Xi|)]/(x 2 ,x 3 ), 

where the first part describes the light particle moving 
around two static heavy particles, and /(x 2 , x 3 ) describes 
for two heavy particles afterwards. By imposing Bcthc- 
Peierls boundary conditions one can find 7 = ±1, and 
the energy of the first part just depends on |x 2 — x 3 |. 
Therefore the wavefunction is reduced to 

^x 1) x a) x a ) = [ V (|r_|)± V (|r + |)]/ 1 (x a -xa)/a(R) ) 

( 51 ) 

and then angular momentum I is determined only by 
/i(x 2 — x 3 ). For 7 = 1, the Fermi statistics require I 
be odd; for 7 = —1, I is even but in this case one can 
easily check that the resultant wavefunction automati- 
cally get rid of the interaction, i.e., Ui/j — 0. Here Fermi 
statistics of two t spins take the crucial role. 

Note that the results presented in (i,ii) uniquely benefit 
from the concept of renormalization and the procedure in 
momentum space to eliminate the ultraviolet divergence. 
These analyses of scattering properties for different mass 
ratios and different angular momenta will be helpful to 
understand the ground state level crossing in the follow- 
ing section. 



2. Energy level crossing 

First, we identify the energy level crossing between 
different angular momenta states for a non-rotating sys- 
tem. The energy spectrum of a non-rotating system was 
previously studied for equal mass[3l|, [32j], and unequal 
masses using Gaussian expansion technique[38j and adia- 
batic hyperspherical method[39]. In Fig^a) we show the 
spectrum for angular momenta I — 0, 1 and for different 
mass ratios using T-matrix method. We also checked for 
higher I > 2 and confirm those states are less modified 
by the interaction and thus not shown here. 
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FIG. 4: (Color online) (a). Energy of three fermions(tti) vs 
interaction strength in a non-rotating isotropic harmonic trap, 
with ao = yl/ijMjf) the confinement length. Different mass 
ratios u — m^/mj. for total angular momentum I = 0(red) 
and / = 1 (black) are plotted. The energy is shifted by the 
ground-state atom-dimer energy, Ed + 3u}/2(Ed is the dimer 
energy). (b)Phase diagram in a rotating harmonic trap with 
rotating frequency SI. The ground state is \l = m = 1) above 
the curves and \l = m = 0) below. 



The system in weak interacting limit(a s — > CP) be- 
haves as non-interacting while in molecule limit (a s —¥ 
+ ) as a single dimer plus an atom. This directly results 
in the inversion of ground state from angular momentum 
I = 1 to I = as l/a s increases. As shown in FigHJa), 
the inversion is denoted by the energy level crossing, and 
the position of level crossing closely depends on mass ra- 
tio u. As expected, when u increases from to 00, all 
energy levels with even-/ move upwards and the system 
evolves from decoupled atom-dimer (except for I = 0) to 
three atoms that are immune from interactions; while 
all odd-/ move downwards until Efimov physics show up 
and invalidate the present T-matrix method. Therefore 
by increasing u, the position of level crossing will move to 
strong coupling side as shown in FigfJIa) ■ Intuitively, one 
can also attribute this to the enhanced s-wave repulsion 
between atom and dimer Q. 

In unitary limit, our numerical results are in good 
accordance with those obtained by using hyperspheri- 
cal coordinates. Previously, hyperspherical coordinate 
method has been applied to a trapped system with equal 
mass[H, |34j. In Appendix IE1 we extend this method to 
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arbitrary mass ratios. Note that the maximum mass ra- 
tio considered in Fig. 3 is much less than critical value, 
= 13.60, for the emergence of Efimov state in I = 1 
channel (as also predicted by Eq |E7l when setting s = 0). 
In this regime, as the matrix size increases we get con- 
vergent result for the energy spectrum. 

Finally, we present the ground state for the trapped 
system with rotation(J7 > 0). By comparing energies 
of all different angular momenta I we obtain the ground 
state as shown in FigfJJb). We find I = m = 1 state is 
gradually favored by the rotation. The energy gain of 
this state is analyzed to be partly from the reduction of 
kinetic energy with respect to I — 0, and partly from the 
avoided s-wave repulsion between atom and dimcr. As 
fi increases, the system evolves to the atom-dimer quan- 
tum Hall state; particularly at O = o>, all states with 
odd-! degenerate. Finally we expect above quantum Hall 
physics of fermionic system could be studied in experi- 
ment, as recently realized in a rotating few-body bosonic 
system (4l|. 



V. SUMMARY 

In conclusion, we present a systematic T-matrix ap- 
proach to solve few-body problems with contact interac- 
tions in the field of ultracold atoms. Taking advantage of 
zero-range interactions, the key ingredient of the present 
T-matrix method is to project the problem to a sub- 
space that is expanded by orthogonal molecular states, 
and meanwhile take careful considerations of the renor- 
malization for relative motions. This method successfully 
unifies the calculations of various physical quantities in 
a single framework, including the bound state solutions, 
effective scattering lengths and reduced interactions in 
the lower dimension. 

We present two applications of this approach, namely, 
two-body scattering resonances in 2D-3D mixed dimen- 
sions and properties of three fermions(tti) m a 3D (ro- 
tating) trap. For the two-body problem, we show that 
T-matrix provides a physically transparent way to under- 
stand the mechanism of induced scattering resonances. 
Besides, it also gives explicit expressions for the reso- 
nance positions and widths. Due to the separate treat- 
ment of relative motions from CM motions, each reso- 
nance can be determined accurately by only consider- 
ing a few related matrix elements. For the three-body 
problem, T-matrix enables us to identify exotic scatter- 
ing properties of three fermions in different angular mo- 
mentum channels and with different mass ratios. In a ro- 
tating system, these properties provide important hints 
for the quantum Hall transition from zero to finite angu- 
lar momentum state. Overall, the external confinements, 
mass ratios, and bosonic/fermionic statistics all play im- 
portant roles and give rise to very rich phenomenon in 
these few-body systems. 
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Appendix A: Construction of individual molecular 
state 

In this appendix, we show how to construct an indi- 
vidual molecular state in the most efficient way. Let us 
consider a system of N particles with masses mi, ...rnjv 
and coordinates Xi, ...x^. For a general case, the molec- 
ular state is written as 

|xi = X2, A = {N, 71,3, ?"i4, ...npf}). (Al) 

In coordinate space it can be factorizcd as 
$at(X) rjil=3 4>ni (xi), where ni (xi) is the eigen-state of 
single-particle Hamiltonian 

V 2 

^o(x,) = -^ + U T ,»(x 4 ) (A2) 
and $tv(X) the eigen-state of 

H (X) = - 7 r r -^ r + Ur. 2 (X) + U T , 2 (X). (A3) 

/(mi + m 2 ) 

Here Vr,i is the trapping potential for the i— th particle. 
The overlap between the molecular state(Eq lAlj) and N- 

particle state(J7j = i is 

N 

(xi = X. 2 ,\\h,l 2 ,h,— In) = lN;h,l 2 JJ^n^j (A4) 

i=3 

with 

1N-MM = J dX$^(X)^ x (X)^(X). (A5) 

The Green function term in Eq llll can then be computed 
efficiently as 

(rj = 0, X\G {E)\rp = 0, A') 
v (rj =0,A|/ 1 ...Z JV )(i 1 .../ jV |r / , =0,A') 
^ E -(E h +... + E lN )+i6 ( > 

ll...(jV 

To this end we have shown a general way to construct 
an individual molecular state. Furthermore, for special 
trapping potentials which enable the decoupling of CM 
motion from other motions, it is convenient to remove the 
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CM motion from the problem and transform the effective 
coordinate vector 



X = (\/2miXi, -v/2m 2 X2, \/2m 3 x 3 , ...^/2m N x N ) (A7) 
to the Jacobi coordinates 



Y = (y/2M^R, y^r, V^IPi, ...^2v N - 2 p N _ 2 ) (A8) 
by a matrix equation 

Y T = AX T , (A9) 
with A-matrix element 



y /rn,j/M N , [i = 1) 

A = J y/A Jj-i/Mj, (£ = j_> l) 



v /m i m J -/(M i _iM i ), (i > j > 1) 
0, all else 



(A10) 



here Mj = £ - = i '• 



For the CM motion, we have the coordinate and the 
mass 



i=i i=i 
for other motions, we take the unique choice as 



2/ir 
2v jP . 



2m\mi 
mi + m 2 



(x 2 - Xl), 



(All) 



(A12) 



/ 2m l+2 M i+ 



V M «+2 

According to Eq lA10| A-matrix satisfies 

AA T = = /, 
here / is identity matrix; this gives 



WE^XAIS) 

i = 1 J 1 



JV-2 



iV 



and 



e?R<ir dpj = J^J dxi 



V X V^ = V Y V^, XX' = YY 



(A14) 



(A15) 



(A16) 



Therefore Mr, /x, v\, ...v^-2 can be considered as 
the effective mass respectively for the motion of 
R, r, p l5 ...p N _ 2 . By decomposing the trapping poten- 
tial to be 

]T VrAxi) = V T (Tl) + Vt(t, Pl , ...p N _ 2 ), (A17) 



we choose the molecular state \r = 0, A) such that its real- 
space wavefunction <j)\{p ll ...p N _ 2 ) is the eigen-state of 
the following Hamiltonian for TV — 2 particles 



H, 



N-2 



N-2 y2 

J2^ + Vr (r = Q,Pi,-P N -2)- (A18) 



The Green function in Eq|TT] is obtained by inserting 
eigen-states of the following Hamiltonian for N — 1 par- 
ticles 



H, 



JV-1 



V 2 N-2TJ-Z 



Compared with Eq. (IAHIA"6l) for a general case, the degree 
of freedom here is further reduced by one particle. 



Appendix B: Derivation of Eq JIOl 

To facilitate the derivation of EqfTUl we assume there 
is only one pair of particles interacting with Ui5 3 (ri), 
and we have 

(r I ,\\U = U I 6 a (r I )(r I = 0,A|. (Bl) 
Using the Lippmann-Schwinger equation or equivalently 
T = U + UG T, (B2) 
and together with EqJ5]we obtain 

£ f y \rv = 0, A') = Ufa) + UG £ fy \n> - 0, A'()B3) 

I'X' I'X' 

The inner product with (n = 0, A/| gives 

T>. 1 



E fx'ijr 6 "' 5 ^' - ^ = °> A l 6 °l r '' = °> A ')) 



I'X' 



(77 = 0,A|V>o>- 



(B4) 



Here we have extracted the most singular terms as char- 
acterized by S 3 (ti). After further renormalizing the bare 
interaction Ui we obtain EqflOl 

Above derivation can be generalized to the case when 
the molecular state \rj = 0, A/) is the superposition of 
many individual ones according to bosonic or fermionic 
statistics. After a proper combination of the resulted 
individual equations, one can equally obtain EqfTOl 



Appendix C: Effective scattering in 2D-3D mixture 

In this appendix we study two-body scattering in 2D- 
3D mixed dimensions(see also Section IIIB). 

First, we relate the effective scattering length to T- 
matrix element between zero-energy scattering states. 
From Eq[TJ the two-body wavefunction in coordinate 
space reads 

^{Pab^AiZb) = iPo{Pab, z a,zb) + 2J zJ <?W0u,ao) 



riA k„,fc 2 



k p ,n A ,k z \T\%), (CI) 



2fi 2rri£ 
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with p AB the relative coordinate of A and B in xy plane 
and z A , z B their respective coordinate in z-direction; E 
is counted from the zero-point energy u> A /2. For large 
separations, both e lkp ' PAS and e lk * ZB oscillate far more 
rapidly in k-space than T-matrix term. Thus in the limits 
of E — > + and |p^b|, \z b \ — > oo, we can specify k p = 
0, k z = in all T-matrix elements and reduce Eq ICll to 



4>(Pab> z a, zb) -> (/)o(z A )(l 



n A >0 



2nd AB 



2nd AB 

-K n a d,A 



A d ^V(0,n A ,0\T\4> ) 



(f) {z A ,a )(l 



Oeff 
d AB 



), (d A B oo) 



(C2) 



with K nA = y/2m B n A uj A , d AB = J^P A B+ Z %' and 
a e ff directly related to T-matrix element as given by 

In Eq|36l e„ and Wn can be obtained from the diago- 
nalization of the following matrix 



Cnn 1 



2irao r 1 



\-Y — 



/(2f) 



kz 



with 



n A kzM p ' E - n A u) A -%L-£i-+iS 



fN;n A ,k z = I dZ<j>* N {Z,a)<f> nA (Z,a )e ik * z '. (C4) 

«/ — OO 

Using the exact identity 

El/^,U 2 = l> (C5) 



2»! 



r, a = y/l/(fiw); L n 2 is 



where Af nl = y r(n+ , +3/2) , 
the generalized Laguerre polynomial; the corresponding 
eigen-energy is E n i m = (2n + I + |)w — mfi. 

The eigen-solution(i?) of the interacting system is de- 
termined by Eq |121 or equivalently (for given {lm}) 



Det(—5 nn > 
a s 



C nn ') — 0, 



(D2) 



where C nn < = A n S nn ' ~ F nn >, with 
f4VA 



lim 

A— >oo 



A 



2 (2n" + l)!! 



n ii—Q E ~ E n lm — E n tiQ0 



F„ 



B,, 



7! E — E n im — E n "Q0 



(D3) 



(D4) 



Br, 



2 f°° 

= (-/3) l Af„iN n n-= / dxx 
V 71 " Jo 



21+2 -x z 



A+\ 



LZ>(x')L^((3V)L*„(a'x'). (D5) 

Here E — E/uj. In our numerical simulations, the typical 
size of the matrix for diagonalization is 50 x 50. When 
sum over intermediate states, the actual cutoff of n", 
namely n c , depends on the values of {n,n'}, which is 
constrained by B nn / n a < 10 in practical calculations. 



Appendix E: Unitary three-fermions with arbitrary 
mass ratios 



we eliminate the logarithmic divergence when summing 
over k p and finally simplify Eq |C3l to be 



Cnn' 



1 f°° f 

- dk z {- \n(k 2 z )S NN , + ^ 
^ Jo „ . 



fN;n Al k z f*N>- nA kz m ( 



1+U 



n A 



^)}C6) 
1 + u J 



In practical calculations, we set the cutoff of n A to be 
1000 for K-Rb{m A < m B ) and 700 for K-U(m A > m B ) 
mixture, depending on the convergence in terms of n A 
in each case. The cutoff of k z (= k z ao) is chosen within 
20 ~ 30 to check the insensitive dependence on the cutoff 
and ensure the accuracy of the integration. 



Appendix D: Evaluation of Eq J45l 

The eigen-state of Hamiltonian (Eq l43|) is 

2 

^nlm(r, 6, </>) = Nnl% m {-) l e- r2 /^L l rt* (^)Y lm (9, , 



(Dl) 



In unitary limit, it is convenient to use hyperspherical 
coordinates to express the three-body wavefunction[Tol 
[H, H3 as 

^(r ± ,p ± ) = ffl)(l-P 23 )£|^r im (p_), (El) 

with r±, p ± given by Eq|40j P23 is the permutation 
operator for identical fermions 2 and 3; 



/rl+pl |r 
R = \l , £ = arctan ±— 



(E2) 



are respectively the hyperradius and hyperangle. Using 
this ansatz, we obtain two decoupled Schrodinger equa- 
tions for R and £, 



h 2 , d 2 1 d s h 2 s 2 



An y dR 2 RdR 4^R 2 



-fiuj 2 R 2 ]F{R) = EF(R): 
(E3) 



d 2 1(1 + 1) 
d£ 2 cos 2 ^ 



W(0 = sM0- (E4) 
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On the other hand, the Bethe-Peierls boundary condition 
applying to Eq lEll gives 



(E5) 



f-iV 

ip'(0) - ^--^(arccos/?) = 0, 

ap 



with a, j3 given by Eq|47l Combined with Eq|E4l EalE5l 
generates the following equations for I = and I = 1 
(assuming all s are real and positive), 

7T 1 

scos(— s) H sin(s arcsin/3) =0, (I = 0) (E6) 

2 ap 



(1 - s 2 ) sin(-s) 



a/3 1 



scos sarcsm 



/3) — — sin(s arcsin^)] 



(1 = 1). 



(E7) 



which are the most relevant equations to the phase tran- 
sition discussed in Section IVB. For equal mass /i = m/2, 
above equations reduce to those in Ref . [HI, HH . 

Note that the validity of above assumption (s > 0) in 
deriving Eqs. (|E6|E7|) closely depends on the mass ratio 
u and the angular momentum I. For any given /, the 



hyperangular equation (IE4I) and condition (IE5[) provide 
the solutions for s, which has double degenaracy(s and 
— s). For I = 1 channel, the zero-range model would 
provide universal energy solutions for u < 8.62 (with 
s > l);42j. When 8.62 < u < 13.6 (0 < s < 1), the 
energy solution will depend on the details of interacting 
potentials and 3-body resonance might happen by tuning 
the potentials [H, Hl| • When u > 13.6 (s 2 < 0) the imag- 
inary s indicates the Ehmov physics with infinite number 
of shallow trimersQ. 



In the actual computation of energy spectrum using T- 
matrix, for nearly all mass ratios below 13.6 one would 
get the energy spectrum with good convergence, as the 
matrix size of Eq |D2l is increased. The signal of 3-body 
resonance(occurs around u = 12.3 with s = 1/2) is van- 
ishing weak due to the infinitesimal width produced by 
zero-range modelpB} . Therefore in the regime of u < 13.6 
one can just consider the s > solution. Then the hyper- 
radius equation (|E3|) gives the energy E = (2n + s + l)u 
with n(= 0, 1...) a semi-positive integer. In the rotating 
frame, the energy is further shifted by —mil, with m the 
magnetic quantum number and fl the rotating frequency. 
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